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NH3 Observations of the Infrared Dark Cloud G28.34+0.06 
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ABSTRACT 



We present observations of the NH3 (J,K) = (1,1) and (2,2) inversion transi- 

> ■ 

O ' tions toward the infrared dark cloud G28.34+0.06, using the Very Large Array. 



Strong NH3 emission is found to coincide well with the infrared absorption feature 
in this cloud. The northern region of G28.34+0.06 is dominated by a compact 
clump (P2) with a high rotation temperature (29 K), large line width (4.3 km 
p |. s~^), and is associated with strong water maser (240 Jy) and a 24 /im point 

O ' source with far IR luminosity of 10^ Lq. We infer that P2 has embedded massive 

c/3 I protostars although it lies in the 8 /im absorption region. The southern region 

has filamentary structures. The rotation temperature in the southern region 
decreases with the increase of the integrated NH3 intensity, which indicates an 
0> ! absence of strong internal heating in these clumps. In addition, the compact 

Q\ . core PI in the south has small line width (1.2 km s~^) surrounded by extended 

emission with larger line width (1.8 km s~^), which suggests a dissipation of tur- 
bulence in the dense part of the cloud. Thus, we suggest that PI is at a much 
O ■ earlier evolutionary stage than P2, possibly at a stage that begins to form a 



cluster with massive stars. 

Subject headings: clouds - ISM: kinematics and dynamics -stars: formation 



^Astronomy Department and CAS-PKU Joint Beijing Astrophysics Center, Peking University, Beijing 
100871, P.R.China 

^Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138, U.S.A. 

■^Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 



-2- 



Introduction 



Stars form from the collapse of dense molecular gas and dust. While low-mass stars 
tend to form in isolation and small groups, high-mass stars ( M > 8 M0) are born almost 
exclusively in clusters (Lada & Lada 2003). The process of low- mass star formation in the 
isolated environment is relatively well studied (Shu, Adams & Lizano 1987). However, the 
processes through which gas and dust condense to assemble massive stars and clusters are 
poorly understood. Klessen (2004) and Mac Low (2004) proposed that, in large scales, the 
supersonic turbulence provides the support to the cloud, while, in small scales, density en- 
hancements cause gravitational collapse. And Myers (1998) and Goodman et al. (1998) also 
proposed that, after the turbulence dissipation, gas condensation will loss internal support 
against the gravity to form stars. Such scenario requires observational tests in massive re- 
gions where clouds are very likely to form massive stars and clusters. Studies of massive 
infrared dark clouds (IRDCs) can provide such opportunities. 

Infrared dark clouds are the dark patches in the sky revealed in the infrared wavelengths 
against the bright galactic background. Surveys along the galactic plane using the Infrared 
Space Observatory (ISO) and the Midcourse Space Experiment (MSX) identified thousands 
of extinction features at 8 /xm (Egan et al. 1998; Hennebelle et al. 2001; Simon et al. 2006). 
Their spatial coincidence with emission from molecular lines and dust indicates that these 
dark clouds consist of dense molecular gas and dust of n^^j ~ 10^ cm^^ (Carey et al. 1998, 
2000). Temperatures in these clouds are low enough (10-20 K) that they do not radiate 
significantly even at the mid or far infrared wavelength. The large amount of dense gas and 
dust in these clouds make them the natural birth place for clusters (Carey et al. 1998, 2000; 
Johnstone et al. 2003; Rathborne et al. 2005, 2006, 2007; Teyssier, Hennebelle & Perault 
2002; Redman et al. 2003; Pillai et al. 2006a,b; Sridharan et al. 2005; Beuther et al. 2005). 

However, most of the previous studies focus on IRDCs that already have strong star for- 
mation activities (e.g. Rathborne et al. 2005), and thus, are at a more evolved evolutionary 
stage. In addition, due to the large distances to these objects, observations with single dish 
telescopes do not have sufficient spatial resolution to reveal the substructures in massive 
clumps and variations of physical parameters. Using the Very Large Array, we imaged a 
giant infrared dark cloud, G28. 34-1-0. 06 at high angular resolution. The object of interest is 
a good target to search for molecular clumps in an early evolutionary stage. At a distance 
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of ~4.8 kiloparsecs (kpc), it contains about 10^ M© of dense gas in the infrared absorption 
region (Carey et al. 2000; Rathborne et al. 2006). The NH3 observations with the Effelsberg 
100m telescope show that, on average, its hne width is smaller than those obtained in high 
mass protostellar objects and Ultracompact HII (UCHII) regions, and the gas temperature 
is lower than 20 K (Pillai et al. 2006b). In our VLA observations, the temperature structure 
and line width variations indicate that the northern region of G28.34 harbors massive proto- 
stellar objects, while, in the southern region, the clumps have characteristics of the earliest 
phase of star formation. 



2. Observations 

We observed the NH3 (J K) = (1,1) line at 23.694 GHz and the (2,2) line at 23.723 GHz 
with the VLA of the NRACU on 2004 June 15 in its D configuration. Seven pointings were 
observed to cover the extended emission; each pointing has an on-source integration time 
of 12 minutes and covers a FWHM primary beam of 2'. We used the 4IF mode that splits 
the 256-channel correlator into four sections to allow simultaneous observations of the NH3 
(1,1) and (2,2) lines with two circular polarizations for each line. The channel separation 
used was 48.8 KHz (~0.6 km s~^ at the line frequency). The time variation of the gains was 
calibrated by quasar, 1851+005, observed at a cycle of ~20 mins. The absolute flux density 
is established by bootstrapping to 3C286. The bandpass is calibrated via observations of 
3C84. 

The visibility data sets were calibrated and imaged using the AIPS software package of 
the NRAO. The average rms of the final images is ~ 3 mJy/beam per 48.8 KHz wide channel. 
The synthesis beam is about 3"x5" when using the natural weighting of the visibilities. In 
order to recover extended structures missing in the interferometer data, we combined the 
visibility data from the VLA with the Effelsberg 100m single dish data from Pillai et al. 
(2006b) for both the NH3 (1,1) and (2,2) lines following the Miriad procedure outlined in 
Vogel et al. (1984). 



^The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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3. Results and Discussions 

3.1. General Properties 

Fig. 1 presents the integrated intensity of tlie NH3 (1,1) line from the combined data set 
overlaid on the 8 /xm image obtained from the Galactic Legacy Infrared Mid-Plane Survey 
Extraordinaire (GLIMPSE; PI: Churchwell) using the Infrared Array Camera (IRAC). The 
dense gas traced by NH3 follows the 8 /im absorption feature very well. Several compact 
clumps are revealed in the filamentary structure. These clumps also coincide with the 1.2 
mm dust continuum emission obtained from the IRAM 30m telescope (Rathborne et al. 
2006). 

In the region north of (5(2ooo) < -04:00:45 (hereafter the Northern Region), both the NH3 
and 1.2 mm emission are centrally peaked with a dust mass of >600 M© (Carey et al. 2000). 
There exists an HIl region associated with IRAS 18402-0403, but has little NH3 or 1.2 mm 
emission. The 24 ^m image from the 24 and 70 Micron Survey of the Inner Galactic Disk 
with MIPS (MIPSGAL; PI: Carey) reveals five embedded point sources. IRS 3 coincides 
with IRAS 18402-0403, and has a luminosity of 10^ Lq (Carey et al. 2000). IRS 2 is detected 
toward the 1.2mm continuum peak P2, and is associated with the NII3 peak and a strong 
H2O maser (240 Jy). Its luminosity of 10^ Lq based on the Spitzer data from 3.6 to 70 /im 
and 1.2mm flux, and a lack of cm emission suggest a high- mass protostellar object. The 
remaining three 24 ^m sources have luminosities of 10^ Lq. 

The region south of (5(2ooo) < -04:00:45 (hereafter the Southern Region) has no detectable 
emission at 1.3cm at an rms of 0.7 mJy, and is associated with three H2O masers (Wang et 
al. 2006). Both the NH3 and 1.2mm emission have a flatter spatial distribution, with a total 
mass over 1500 M© (Rathborne et al. 2006). Four 24 /im sources are detected from the MIPS- 
GAL survey with luminosities < 10^ Lq. The strongest source, IRS 1 (q;(2ooo) = 18-42:50.82 
^(2000)= -04:03:09) with a luminosity of 10^ Lq, is associated with an NH3 peak toward 
PI (a(2000) = 18:42:50.27 (5(2ooo)= -04:03:20), and an H2O maser emission (a(2ooo) = 18:42:50.80 
'^(2000)= -04:03:11.0). The other two H2O masers in the Southern Region has no 24 /xm 
sources within a projected distance of 30". 

It is difficult to predict whether the Southern Region will eventually form 
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massive stars similar to the Northern Region. The Southern Region contains more 
dense molecular gas than the Northern Region. Assuming a similar star forming efficiency, 
one would expect a more massive cluster to emerge in the Southern Region. However, 
since the south is spatially more enxtended, it may form less massive stars than 
the Northern Region. On the other hand, if the cluster follows a normal Salpeter 
IMF, one would expect more massive stars in the Southern Region. Thus, a lack 
of luminous infrared sources in the south indicates that it is at an earlier evolutionary stage 
than the Northern Region. Furthermore, the radio data in conjunction with the infrared data 
appear to suggest an evolutionary sequence of star formation, starting with an evolved stage 
in IRAS 18402-0403/IRS3; to an embedded phase in IRS2 which has a luminosity similar 
to high-mass protostellar objects; to less luminous IRSl which has strong mm continuum, 
and H2O maser emission; and finally to NH3 clumps which only has strong mm continuum 
emission, but no 24 /im emission peak nor H2O maser emission. 

The NH3 data reveal physical properties that further conffim the difference in the North- 
ern and Southern Region. In the Northern Region, there is a dominant compact clump P2 
associated with IRS 2 (a;(2ooo) = 18-42:51.92 5(2ooo)= -03:59:54), which has a peak rotation 
temperature close to 30 K and a FWHM line width larger than 3 km s~^, and is surrounded 
by an extended envelope. By averaging rotation temperatures in the Northern Region of 
the same NH3 (1,1) intensities, we derived a relation between the two quantities. Fig. 2 
presents the average NH3 temperature within an intensity bin versus the mean value of the 
bin. As shown in the black line in Fig. 2, the rotation temperature in the Northern Region 
increases with the increase of the integrated intensity, which indicates strong internal heating 
in clump P2 of the Northern Region. Therefore, although clump P2 is embedded in the 8 yum 
dark region, the high temperature, large line width, strong water maser, and 24 /im point 
source indicate that this compact clump is forming massive protostellar object (s), similar to 
the compact cores found in other IRDCs (Pillai et al. 2006a; Rathborne et al. 2005, 2007; 
Beuther et al. 2005). 

On the other hand, the NH3 emission from the Southern Region appears to have different 
characteristics as compared to that in the Northern Region. The red line in Fig. 2 shows that 
the rotation temperature in the Southern Region decreases with the increase of the integrated 
NH3 intensity. On average, areas with stronger NH3 emission have lower temperatures, which 
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indicates absence of strong internal heating in the southern clumps, unlike P2 in the north. 
Toward the most compact clump Pl/IRSl in the Southern Region, the NH3 gas has a lower 
rotation temperature (~16 K) and narrower line width (~1.8 km s^^) than that in P2. We 
also note that the temperature decrease appears to flatten toward the highest intensity in 
the red line in Fig. 2. This could be due to the limited heating from star formation activities 
in the area. 

We also compared the general properties of G28 Southern Region with the local, in- 
termediate mass star formation region OMC-north (Bally et al. 1987; Hillenbrand 1997; 
Castcts & Langer 1995; Johnstone & Bally 1999). The linear resolution of our interferome- 
ter and single dish combined NH3 image (~4" resolution) of G28.34+0.06 (distance of ~4.8 
kpc) is comparable with the Effelsberg 100m single dish NH3 image (40" resolution) of the 
OMC-north (distance of ~0.48 kpc) (Cesaroni et al. 1994). We found that, on average, 
G28.34 Southern Region have lower temperatures (16 K vs. 20 K), larger line width (1.6 vs 
0.8 km s^^) and higher column density (2.3 vs. 0.78x10^^ cm^^) assuming the same NH3 
abundance (0.3 x 10^^) in both regions. The mass within the dense filament estimated from 
dust emission is ~1500 M0 in G28.34-south region (Rathborne et al. 2006), while the dense 
mass in OMC-north is ~350 M© (Lis et al. 1998). On the other hand, G28.34-south region 
extends 4' (~6 pc) at the distance of ~4.8 kpc, while OMC-north extends 20'(~3 pc) at the 
distance of ~0.48 kpc. All these parameters indicate that G28.34-south has colder and more 
turbulent dense gas packed in a smaller region than OMC-north. 



3.2. Line Width Variation of NH3 Emission toward Peaks PI and P2 

We examined the spatial variation of NH3 (1,1) line width from the VLA data as well 
as that combined with the Effelsberg 100m single dish data. The line width variations after 
the data combination provide kinematic information on the scale of > 0.1 pc (4" in 4.8 kpc). 
The line width of PI and P2 before and after the data combination are shown in Table 1. 
It appears that the line width in PI is broader after adding the single dish data, while, in 
P2, we find an opposite effect. Since the missing flux in an interferometer corresponds to 
the extended spatial structures, the variation in the line width represents a change in the 
level of turbulence in their envelopes/ cores along the line of sight. Normally, in a typical 
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active massive star forming core, such as P2, the missing fluxes in an interferometer are 
from extended and quiescent envelopes. Dense and compact emission surrounding high mass 
protostars or UCHII regions has broader line width due to motions such as infall, outflow 
and rotation. The missing extended structures are normally weaker and have narrower line 
width. However, toward the G28.34 Southern Region, we observed the opposite effect: the 
extended emission shows systematically larger line widths than those from the compact core 
component revealed by interferometer. 

This effect is also shown in line width distribution in Fig. 3. The line width distribution 
appears to behave differently in the clumps from the Northern and Southern region after data 
combination. As shown in Fig. 3, for the cold clump PI in the Southern Region, the mean 
line width becomes broader after recovering the missing flux, while for the warm clump P2, 
the mean line width is similar to the one from the interferometer data alone. Furthermore, 
the line widths of the entire Southern Region follow a distribution similar to PI. We flnd 
that the effect of relative motions among clumps can be ruled out. We tested this scenario 
by convolving the combined interferometer data to the resolution of the single dish telescope 
(40"), and comparing with the line width of the combined data at 4" resolution. The mean 
line widths from these two data sets are similar. 

The variations of NH3 line width in the G28 Southern Region suggest a dissipation of tur- 
bulence, possibly at the onset of star/cluster formation. The dissipation of turbulence could 
initiate gravitational collapse as suggested by several observations (e.g. Ward-Thompson, 
D., 2002). Taking advantage of the high resolution of interferometers, we obtained the line 
width variation of PI and P2 along the line of sight on scales of >0.1 pc. Using the NH3 
data, we calculated the molecular mass and virial mass (Estalella et al. 1993) of PI and 
P2 (Table 1). In PI, the material seen from both interferometer alone and the combined 
data sets is gravitationally bound. And the dense component is more likely to collapse (M 
= 16 Mq vs. M^ir = 6 Mq) than the component including all the spatial structures (M = 
27 M0 vs. Mmr = 14 Mq). lu coutrast, due to strong star formation activities, the dense 
component in P2 has a larger virial mass (75 Mq) than the molecular mass (30 Mq), while 
the whole clump is gravitationally bound and may undergo collapse (M = 75 Mq vs. M^ir 
= 44 Mq). We suggest that, in the cold massive clump PI, and the rest of the G28.34-south 
region, dense gas undergoes fragmentation and condensation. Considering the low rotation 



temperature (~16 K), relatively narrow line width (~1.8 km s~^), PI appears to be in a 
very early evolutionary stage during which turbulence dissipates to facilitate gravitational 
collapse to form stars. At this stage of cluster/massive star formation, the most massive 
members in the cluster may be still at the intermediate mass stage, and will likely continue 
to accrete to become more massive stars, given the amount of dense gas available in the 
region. 
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Fig. 1. — The integrated intensity of the combined NH3 (1,1) emission in white sohd contours 
overlaid on the Spitzer 8 //m image in logarithmic color scales. The NH3 image is contoured 
at 10% of the peak (1 Jy beam^^ x km s~^). The star symbols mark the 24 fira emission 
peaks observed with MIPS/ Spitzer. The cross symbols mark water maser emission detected 
with the VLA. The thin dashed line indicates the 50% of the sensitivity level of the 7 pointing 
mosaic in NH3. The NH3 data have a resolution of 5" x 3" , shown as the shaded ellipse at 
the lower-left corner of the panel. 
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Fig. 2. — Relation between rotation temperature and NH3 (1,1) integrated intensity in the 
Northern Region (black) and the Southern Region (red) of G28.34+0.06. The error bars 
present the la standard deviation. 
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Fig. 3. — NH3 (1,1) line width distribution in different velocity bins toward G28.34+0.06. 
Red columns are the number count from the VLA data set. Green columns are the number 
count from the combined data set (VLA+lOOm). Left: Line width distribution for clump 
P2 over an area of 10" x 10". Right: Line width distribution for clump PI over an area 
of 15" X 5". The NH3 line widths, measured in FWHM, are corrected for optical depth and 
broadening due to hyperfine structures. 
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PIC 


18 42 50.27 


-04 03 20 


1.2/1.8 


0.59/1.0 


16/27 


6/14 


P2<= 


18 42 52.07 


-03 59 54 


4.3/3.3 


0.48/1.2 


30/75 


75/44 



Table 1. Physical Parameters of PI and P2 over the Scale of 0.1 pc 

Source R.A. (2000) DEC. (2000) Aii" N^oi" Mass'' Virial Mass" NH3 Abundance'' 

Name (km s^^) (lOi'^cm-^) (Mq) (Mq) (10^^) 

0.6 

0.3 



Note. — Units of right ascension are hours, minutes, and seconds, and units of dechnation are degrees, arcminutes, 
and arcseconds. 

^The first number is obtained from the VLA data only; The second number is obtained from the combined data. 

''The NH3 abundance is calculated by comparing the column density of dust (Rathborne et al. 2006) and NH3 
emission. We convolved the NH3 data to 11", the the resolution of the dust emission, and assumed that the NH3 
abundance is the same in (11" X 11"). 

•^Thc coordinates of PI and P2 refer to the integrated intensity peaks of combined NH3 (1,1) data. 



